Enhancement of mobilities in a pinned multidomain crystal 
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Mobihty properties inside and around degenerate domains of an elastic lattice partially pinned on 
a square array of traps are explored by means of a fully controllable model system of macroscopic 
particles. We focus on the different configurations obtained for filling ratios equal to 1 or 2 when 
the pinning strength is lowered. These theoretically expected but never observed configurations are 
degenerated, which implies the existence of a multidomain crystal. We show that the distinction 
between trapped and untrapped particles that is made in the case of strong pinning is not relevant 
for such a weaker pinning. Indeed, one ought to distinguish between particles inside or around 
the domains associated to positional degeneracies. The possible consequences on the depinning 
dynamics of the lattice are discussed. 

PACS numbers: 74.25. Qt 



I. INTRODUCTION 

When submitted to an external force, an elastic lattice 
of interacting particles coupled to a regular array of pin- 
ning centers can move as soon as the applied force exceeds 
a threshold. The determination of this threshold and of 
the behavior of the elastic lattice near it are important 
and current questions. This threshold essentially depends 
upon two parameters : the local deviations of the parti- 
cles with respect to their positions in the unpinned lat- 
tice, and the local energetic landscape around the pinned 
equilibrium configuration. Thus, knowing how this de- 
formation will take place and exploring how the particles 
locally move around their pinned positions in equilibrium 
configurations is a first step towards the understanding 
of the lattice depinning process. 

In this article, we show how the existence of multido- 
main configurations would act on the dynamics of an 
elastic lattice pinned by a periodic potential of medium 
amplitude. In order to explore this question, we have 
developed a novel experimental system allowing the di- 
rect observation of the interacting particles and a well 
controlled pinning process. It will be also presented. 

Many elastic systems submitted to a pinning array 
have been explored for some years. The most studied 
of them are the flux lattices in superconductors since a 
good pinning of the vortices increases the critical current 
in the material. Recently, the focus on such questions 
has grown since progress in lithography techniques on 
superconducting thin materials has offered the opportu- 
nity to design arrays of holes or magnetic dots, which act 
as strong pinning centers where one vortex or even more 
can be trapped -^^^diL^ this procedure is much more 
efficient and potentially controllable than the weak pin- 
ning resulting from intrinsic defects in the material^i In 
such a situation, a key parameter is the ratio / between 
the vortex density, which is controlled by the magnetic 
field intensity, and the pinning centers density. Indeed, 
measurements of macroscopic characteristics such as the 
voltage-current response or of the vortex lattice melt- 
ing temperature have shown unambiguously that they 



strongly depend upon / for a square or triangular pinning 
array and that the pinning is more efficient for rational 
p/q values of /, where p and q are small ^ 

The dependence of the depinning threshold upon the 
microscopic vortex organization and their dynamics has 
been proved through semi-analytic or numerical works 
that were first devoted to the static problem of the pinned 
lattice configuration j^i-'^^i-'^-'^i-'^^i-'^^i-'^^ Taking the lowest en- 
ergy equilibrium configuration as a starting point, the 
different phases of a driven moving lattice have been sim- 
ulated |11i1^i15t16 ^YiQ resulting quasi- ID movement has 
also been specifically studiedJ^ Let us emphasize at this 
point that, whereas in most cases the equilibrium con- 
figurations are degenerate, the possible existence of mul- 
tidomain lattices is seldom taken into account in these 
numerical studies We shall see that it could be a prob- 
lem since particles' behaviors strongly differ depending 
on whether they are inside the domains or in the walls 
between them. 

From the experimental point of view, investigating 
the microscopic configurations of a vortex lattice is still 
a challenging issue. The best imaging of equilibrium 
configurations has been obtained with scanning Hall 
probe,^-5^i^i^ scanning tunneling^ or Lorentz^ microscopy. 
In these works, configurations of mainly some hundreds 
of vortices are imaged for rational filling ratios / close 
to 1 (1/4 < / < 4). However, even though great efforts 
have been made in this way, a direct imaging of (pinned) 
vortex lattices suffer from two major drawbacks. First, 
it is presently impossible to track the vortices and follow 
their dynamics ; furthermore, the nature of the artifi- 
cial pinning centers used in the works cited above does 
not offer the required experimental facility to control the 
amplitude and the shape of the pinning potential. 

Therefore, a complementary approach based on simi- 
lar systems has appeared to be necessary. Beyond their 
intrinsic interest, colloidal crystals can be seen as an at- 
tractive alternative solution. The first attempts in this 
way have recently involved colloidal crystals pinned by 
arrays of optical traps. Nevertheless, the observed con- 
figurations, which are also numerically predicted^i^S are 
to date either quite different than the one observed in 
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FIG. 1: Summary of the two electrostatic forces at stack, (a) 
: Scheme accounting for the origin of the repulsive interac- 
tion between two balls and picture of the resulting triangular 
lattice, (b) : Scheme accounting for the apparition of the pin- 
ning force under two areas linked to two different potentials 
and picture of a few adjacent patterns in the ITO electrode 
that is used to study the f = 1 case. 
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FIG. 2: Overview of the equilibrium configurations on a 
square array with filling ratios equal to 1,1/2, 1/3 and 1/4. 
The corresponding theoretical schemes for a mono-domain 
configuration are also shown. The traps are represented by 
empty circles, while the black disks represent the particles. 
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superconductors,^^ or restricted to system of very small 
spatial extent. Moreover, the role of the hydrodynamic 
couplings and of the nature of the pinning strength are, 
in these systems, still unclear. 

In order to bypass these difficulties, we propose here 
an alternative experimental system consisting of a few 
thousand millimetric charged particles which can be sub- 
mitted to a controlled electrostatic pinning potential and 
a tunable effective temperature resulting from a mechan- 
ical shacking. After a brief description of the pinning 
set-up and the experiments which validate its efficiency 
(Sec. |TT]), we present in Sec. Illll the multidomain equilib- 
rium configurations and their evolution when the ampli- 
tude of the pinning strength decreases in the case of a 
square array of traps. We mainly focus on the / = 1 and 
/ = 2 cases. Finally, the thermally activated motion of 
the particles in these multidomain systems are described 
and discussed in Sec. [iVl 



II. EXPERIMENTAL SET-UP 

The system we propose, a so-called "pinned macro- 
scopic Wigner crystal", is made of a mono-layer of con- 
fined charged millimetric balls. These conducting balls 
are located on the grounded bottom electrode of an hor- 
izontal plane condenser. When a tunable voltage Vi > 
(Vi 1 kV) is applied between the top and the bottom 



FIG. 3: Overview of the equilibrium configurations and cor- 
responding schemes for a triangular array with filling ratios 
equal to 1,1/2 and 1/3. 



electrodes, the balls get negatively charged and repel each 
other (Fig. [Tfa)). In Ref. 23, we have shown that the in- 
terball interaction potential E{r) is equal to EoKo{r/X), 
where Kq is the modified Bessel function of the second 
kind, and Eq is proportional to . This interaction be- 
tween the balls has the same space-dependence as the 
interaction between vortices in superconductors, which 
suggests that the observed behaviors could be directly 
mapped onto vortex systems without the experimental 
constraints associated to the real ones. To introduce ther- 
mal noise, the whole cell is fixed on loudspeakers supplied 
by a white noise voltage. We have thoroughly checked 
that the resulting shaking of the balls, due to friction with 
the bottom electrode, fulfil the properties one can expect 
for a thermal shaking. First, even though all the balls lie 
on the same solid substrate, their movement is spatially 
non correlated, which might be due to inhomogeneities 
on the wafer at a microscopic level. Secondly, the indi- 
vidual trajectory of a single ball which is free or trapped 
in a parabolic well can be described through Langevin 
formalism. Stationary states are reached in a few tenth 
of second. They are characterized by an effective temper- 
ature directly controlled by the shaking amplitude : the 
energy distribution is given by Boltzmann statistics. This 
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was proved on confined small islands of balls, that can 
be seen as two-level systems when considering their two 
first equilibrium configurations, characterized by concen- 
tric shells of varying number of balls. This effective 
temperature was calibrated and is measured in-situ. 

When submitted to a low effective temperature, the en- 
semble of balls can organize and form a defect-free trian- 
gular lattice of a few thousand particles By insulating a 
periodic set of small patterns in the top electrode we can 
link it to a different tunable voltage V2 > Vi, and create 
an array of pinning sites, as the local potential difference 
results in an attractive force on a ball (Fig. [TJb)). Since 
the net charge on the ball around the trap is roughly pro- 
portional to the mean potential V between Vi and V2 , 
the amplitude of the pinning potential can be considered 
as proportional to VAV^ where AV = V2 — V1 comes from 
the gradient of the electrostatic potential between inside 
and outside the trap. Therefore, the relevant parameter 
to discuss the relative strength of the pinning compared 
with the elastic stiffness wih be r]v = VAV/V'^ = AV/V, 
which is fully controlled. 

In our experiments, the lattice spacing ao in the un- 
pinned lattice (made of nearly 2000 balls) is 1.82 mm, 
to be compared with the screening length A = 0.48 mm. 
The periodicity of the square pinning array is then cho- 
sen so that the filling ratio / is equal to 1 or 2. The width 
of the patterns is such that the range of the pinning po- 
tential associated to one site is of the same order as the 
interball distance. 

In order to test first the efficiency of this pinning device 
we have observed the equilibrium configurations obtained 
when starting from / = 1 and taking off some balls so 
that / < 1, with high values for the relative pinning 
strength {rjv > 0.25). These configurations have been 
compared to those reported in the literature. 

Our experimental equilibrium configurations are 
shown for some fractional filling ratios on Fig. [2] for a 
square pinning array and on Fig. [3] for a triangular array. 
In the first case, they are similar to those observed in su- 
perconductors^^^^^ and numerically predicted. ^^^^^ In the 
second case, they are also similar to those numerically 
predicted. These results confirm the ability of our sys- 
tem to be generic of elastic systems and to describe well 
the behavior of flux lattices in superconductors. 



III. MULTIDOMAIN EQUILIBRIUM 
CONFIGURATIONS 

The existence of multidomain equilibrium configura- 
tions depends upon different parameters. 

For / < 1, the equilibrium configurations are struc- 
turally degenerate since the organized pinned lattice can 
be pinned on various subarrays of traps whatever the 
shape and the intensity of the pinning potential. There- 
fore it is experimentally difficult to obtained large coher- 
ent domains and multidomain systems are often observed 
as for instance in our system (Fig.[4|), or in superconduc- 
tors, where a multidomain experimental configuration of 




FIG. 4: Snapshot of the crystal in a quasi-equilibrium state 
with / = 1/2 and a square pinning array. The theoretical 
configuration is 2-fold degenerate, which leads to the existence 
of two kinds of domains (black or gray background colors). 



a / = 1/2 crystal of more than 10^ vortices has clearly 
been imaged.^ On the contrary, according to what stated 
in Ref. 11, it seems to be difficult to find these multido- 
main configurations through numerical simulations. 

Beyond these easily predictable cases, theoretical stud- 
ies suggest that for / > 1, lowering the amplitude of the 
pinning strength or widening the traps can lead to par- 
tially pinned (PP) configurations, where only a fraction 
of the vortices are located in the middle of the pinning 
sites, allowing multidomain configurations. Such con- 
figurations have never been observed in superconductors, 
perhaps because the pinning intensity is too high ; how- 
ever our experiments in which the pinning strength can 
be tuned confirm this assumption. 

For / = 1, a sharp transition between the totally 
pinned (TP) and the partially pinned configuration is 
observed when lowering rjv ij] < 0.15). The latter con- 
figuration consists in an alternation of rows of trapped 
or untrapped balls, so that the triangular symmetry of 
the unpinned lattice is almost restored, the position of 
the untrapped particle between two free traps slightly 
depending on the relative pinning strength (Fig. [5]). Let 
us now emphasize that this partially pinned configura- 
tions is degenerate since there are two possibilities for 
the orientation of the unpinned rows (Fig. [6]). 

Note that we don't discuss here the existence of a 
"floating crystal", a triangular lattice slightly distorted 
by the traps : in the literature, its orientation is contro- 
versial, as it can be aligned with the axes of the underly- 
ing pinning array or slightly rotated Experimentally, 
this orientation should depend strongly on the boundary 
conditions, unlike the pinned configurations. Anyhow, as 
for / = 1,^^ one should expect this configuration to be 
favored when the traps widen and can include many par- 
ticles : then, a triangular-like structure will be obviously 
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FIG. 5: Thermally activated trajectories around equilibrium configurations and corresponding schemes. For / = 1, the transi- 
tion between both configurations is obtained by tuning r]v- The Delaunay triangulation of the partially pinned configuration 
for / = 2 is also shown. 



preferred. 

In the / = 2 "totally pinned" configuration 

(TP), where all the traps are occupied and the remain- 
ing particles are located in the center regions between 
the traps is usually observed in superconductors^^ and 
also numerically predicted.^ We have also obtained this 
configuration with a strong pinning (77 > 4). For traps 
with larger extent totally pinned configurations switch 
to a regular array of dipoles centered on the traps as 
it was shown by numerical simulations for vortices with 
logarithmic interaction.^^ This configuration could be the 
one observed by Field et al., where two flux quanta are 
trapped in each pinning center."^ Similar configurations 
seem to be also observed in colloidal systems. Indeed, 
the latter affirmation is an extrapolation of the observa- 
tions made in Ref. HJ, where the / = 3 case is studied 
and an array of trimers is observed. This extrapolation is 
confirmed by numerical simulations made by Reichhardt 
et al.^^ In these cases, the degeneracy is linked to the 
rotation of the dimers, while the TP configuration is not 
degenerate. These results indicate that the existence of 
multidomain systems depends upon the intensity and the 
shape of the pinning potential. 

Since in our experiment we can tune the pinning in- 
tensity (the trap shape being fixed) we have explored the 
evolution of the TP configuration as this intensity de- 
creases. For not too high values of T^y, we have observed 
an unexpected and interesting partially pinned (PP) con- 
figuration shown on Fig. [5l Only half the traps are oc- 
cupied while the other ones are the centers of trimers, 
effectively trapping three particles instead of one. It re- 
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FIG. 6: Thermally activated trajectories in a multidomain 
partially pinned configuration for / = 1. 



suits from this an elegant four-fold configuration in which 
each particles has got exactly five equidistant neighbors. 
These patterns result from a compromise between the 
symmetry imposed by the pinning array and the sym- 
metry of the unpinned lattice, which is six-fold. Since 
there are two possibilities for the choice of the subarray 
of occupied traps and, for each selected subarray, four 
possibilities for the orientation of the trimers, the system 
is highly degenerated and, as a consequence, a multido- 
main configuration will appear in the pinned lattice. 



5 



Clir^CK^CikyiDiSD^V4i^ii>Aiirj 



0.10 



I 



II 



HI 



FIG. 7: (color online) Thermally activated trajectories in a 
multidomain partially pinned configuration for f — 2. The 
traps are located at the intersections of the red lines. The 
green circles mark one of the two subarrays that can be oc- 
cupied by particles. 
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FIG. 8: Phase diagram for the three pinned configurations 
that are found here and in the literature for f — 2 and a square 
pinning array. The values for ao and A are the experimental 
ones. The profile of the wells is gaussian, with a depth Ep 
and a width a. 



The multidomain partially pinned configurations are 
observed in our whole system as shown on Fig.[6]for / = 1 
and Fig.Ofor / = 2, where three domains are present. In 
the latter, the difference between domains I and II lies in 
the choice of the occupied subarray. Between domains II 
and III, the trimers do not have the same orientation. In 
both cases, the walls are marked by features looking like 
"eyes" constituted by a well trapped particle surrounded 
by very mobile ones. We shall return to this last point 
in the following. 

In these configurations, the domains' sizes increase for 
high V, thus for high intensities of the forces at stack. 
Notice furthermore that these configurations are highly 
symmetric, the walls between the domains being mainly 
orientated along the principle axes of the pinning lattice. 
These equilibrium multidomain configurations are very 
stable and can be considered as quasi-static because of 
the presence of pinned particles. The evolution in shape 
and size of each domain is much slower than in the case 
of a disoriented cluster in an unpinned lattice that was 
discussed in Ref. 126|. 

In order to determine the conditions required to ob- 
serve such configurations, we have calculated the phase 
diagrams of the pinned macroscopic Wigner lattices, the 
two relevant parameter being the relative pinning inten- 
sity Ep/E^ and the width a of each electrostatic trap. 

Since the analytical determination of these attracting 
potentials is a tricky problem, we have calculated the 
phase diagram for different profile and compared them 
to the experimental results. 

This procedure is relevant to find the depth Ep and 
width cr of the traps since the obtained phase diagrams 
strongly differ according to the choice in the profile used 
for its calculation as it was shown in Ref. [3 for / = 1 

Various shapes have been proposed in the literature : 



in the numerical simulations of flux lattices this poten- 
tial is either a truncated parabola or a gaussian, ^-^^^liiii^ 
whereas the optical traps in colloidal systems seem to be 
well described by sine functions. Our calculations show 
that a good agreement is obtained for traps with gaussian 
profiles Eg{r) = -Ep exp [ - r'^/{2a'^)] . 

In particular, using this trap profile, we have calcu- 
lated the phase diagram for / = 2 for which we consider 
three possible configurations : the totally pinned config- 
uration, the partiall y p inned one and the dipole configu- 
ration found in Ref. [lO|. The resulting phase diagram for 
a gaussian profile is shown on Fig. [HI 

It can be seen from the diagram that the partially 
pinned configuration is favored relatively to the totally 
pinned one for sufficiently wide or shallow traps. Per- 
haps the dipole configuration could not be observed in 
our system because the traps are not wide enough, or 
because rjv has a upper bound by construction. 

Let us nevertheless indicate that such a phase diagram 
has to be considered more as a guide for the experiment 
than a tool for precise predictions : indeed these phase 
diagrams are sensitive to the details of the unknown pro- 
file of the traps and one cannot assert that the real traps 
in our experiments or in flux lattices can be described by 
the simple functions that are usually tried. 



IV. PARTICLES DYNAMICS BETWEEN AND 
ALONG THE DOMAINS WALLS 

The presence of quasi-static domains leads to the ex- 
istence of two kinds of particles, viz. those inside the 
domains, and those near the walls. As discussed by Grig- 
orenko et al.,^ lacking or excess particles are required for 
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FIG. 9: Schemes of the main movement of the particles in- 
side a domain or along a domain wall at low temperature. 
Note that the lacking or excess particles of the wall are not 
represented. 




FIG. 10: Trajectories inside an homogeneous domain of the 
partially pinned configuration for f — 2. In (b) the tempera- 
ture is twice the temperature in (a). 



the adjustment of the domains and the construction of 
the walls. Thus, these two kinds of particles are associ- 
ated to different local densities. Our experiments show 
that even if these density variations are tiny, they induce 
important differences in the dynamics of the balls, even 
in the absence of an external force. More precisely, the 
movements in or alongside the domains strongly differ 
from the view of the associated mobility and the direc- 
tion of the movement. 

In the / = 2 case (Fig. [7]), if one first looks inside 
the domains, the trapped or untrapped balls (viz. those 
belonging to a trimer) are stable and their dynamics es- 
sentially result from the thermal activation. As it can be 
seen on the trajectory spots, the amplitude of the move- 
ments of the trapped or of the trimer particles are equal. 
We can conclude from this result that the particle status 
towards the traps is not relevant to discuss its dynamical 
properties. Underline that this behavior can be associ- 
ated to the conditions under which the partially pinned 
configuration is energetically favorable, notably the ne- 
cessity to have wide traps. In already studied totally 
pinned configurations ^i^^ interstitial particles are much 
more mobile than trapped ones. 

The dynamics along the walls is pretty different. In- 
deed, the "eyes" in the walls are constituted by a well 
trapped particle, whose wandering amplitude is similar 
to the one inside a domain, surrounded by much more 



mobile untrapped particles. This high mobility can be 
associated to the local metastability linked with a small 
density variation. This makes the fact that the trapped 
particles remain in their place and are not chased by their 
neighbors all the more surprising. 

Furthermore, the favored directions are also pretty dif- 
ferent according to the class of the particles. These direc- 
tions are schematized on Fig. [9l Inside a domain, all the 
movements have a privileged direction following the sym- 
metry of the underlying pinning array. Notice that for 
the untrapped balls, this movement is surprising since it 
is not orientated towards the traps. This appears clearly 
when the temperature increases, the symmetry of the 
underlying pinning array being then directly recovered 
through the trajectories mapping (Fig. [T0|). Note that 
this is completely different than the "two-stage melting" 
that is found in the dipole configuration.^^ On the other 
hand, the untrapped particles near a metastable area are 
clearly following trajectories linking two adjacent empty 
pinning centers. We think that this phenomenon should 
be directly linked with the existence of trimers which 
are in a very uncomfortable position and are highly per- 
turbed by a density variation. By contrast, for / = 1, 
the movements along the walls are also enhanced rela- 
tively to the movements inside the domains but are still 
respecting the same 4-fold underlying symmetry as inside 
the domains (Fig. [6|). 

As a conclusion, this study performed with a macro- 
scopic Wigner crystal has shown without ambiguity that 
the existence of domains must be taken into account to 
discuss the particles dynamics in a pinned elastic lat- 
tice. Those multidomain configurations appear as soon 
as the pinning amplitude diminishes and implies degener- 
ate equilibrium configurations. The tracking of the par- 
ticles has shown that the discrimination that should be 
made between trapped and interstitial particles in the to- 
tally pinned case is not relevant in a weaker pinning case 
for which one should discuss the dynamics by discrim- 
inating between particles in stable or metastable situa- 
tions, respectively inside or at the boundaries of the do- 
mains. This suggests that the existence of domains must 
be taken into account in further experimental or numeri- 
cal studies of the moving phases of a partially pinned lat- 
tice submitted to a depinning force. Consequences could 
be first an easier depinning along the walls, which would 
then induce anisotropy since the orientations of the walls 
are given by the principle axes of the pinning array. In 
return, depinning of the lattice could remove the degen- 
eracies, which would then imply isotropy. In Ref. 11, the 
nature of the flow is discussed according to the values of 
the filling ratio /. It is shown that if / is rational, this 
flow has elastic properties, whereas plasticity is invoked 
when / is irrational. Even though the "macroscopic fill- 
ing ratio" is rational here, the local ratio near the bound- 
aries is irrational and thus, the nature of the depinning 
dynamics could become very complex with a coexistence 
of elastic and plastic flows in the same lattice. 
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